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Abstract 

Liquid is represented as complicated systera of disclinations ac- 
cording to defect description of liquids and glasses. The expressions 
for the linear disclination field of an arbitrary form and energy of inter- 
disclination interaction are derived in the framework of gauge theory 
of defects. It allows us to describe liquid as a disordered system of 
topological moments and reduce this model to the Edwards- Anderson 
model with large-range interaction. Within the framework of this ap- 
proach vitrifying is represented as a "hierarchical" phase transition. 
The suggested model allows us to explain the process of the Fischer 
clusters formation and the slow dynamics in supercooled liquids close 
to the liquid-glass transition point. 

1 Introduction 

The Fischer clusters are long-range correlations of density fluctuations, which 
are observed in the supercooled liquids (~100 K above the glass transition 
temperature Tg) [1]. To date these clusters have been discovered in the 
polymers, glassforming, as well as in single fluids. The typical size of these 
formations is ~ 100-300 nm, whereas the correlation radius of usual short- 
range thermal density fluctuations is ~1 nm. The fractal- like structure and 
long lifetime are important typical properties of the Fischer clusters [2]. 
There are theoretical works in which the formation of these clusters is sup- 
posed to be the process of condensation of atom groups with a common 
type of local atomic ordering [2, 3]. However, in the case of dense liquids 
(at low temperatures) in these models have complications in definition of 



order parameter. Besides this approach does not ahow us to explain the 
roots of fractal-hke structure of these formations. Therefore the nature of 
the Fischer clusters is not quite clear yet. 

In the present work the theoretical model of glass transition, based on 
the known approaches of the disordered systems physics, is offered. Our 
model is based on the disclination model of amorphous structures, sug- 
gested by D. Nelson. Prom the beginning we write the system Lagrangian 
and deriv the expressions for the linear disclination field and energy of the 
inter-disclination elastic interaction (2-th section). After that we define the 
topological moment of the disclination system and describe liquid as a dis- 
ordered system of interacting topological moments. It allows us to represent 
the liquid-glass transition in terms of spin-glass systems physics, and use 
the spin-glass theory to estimation of the size of the long-range correlations 
in supercooled liquids (3-th section). 

2 Disclination description of liquids 

The disclination description of liquids, suggested by D. Nelson, is the basis 
of our model. The basic thesis of this one is that the regular tetrahedron 
represents the closest and most profitable local atomic packing. Indeed, it 
has been demonstrated that the structure of atomic liquids and glasses has a 
significant polytetrahedral character [4], first due to the success of the close 
random packing of hard spheres [5] as a model for metallic glasses and later 
due to computer simulations [6]. In the framework of this polytetrahedral 
model the liquid structure is considered to be a tessellation of all space by 
tetrahedra with atoms at their vertices. However, it is known, that Euclidean 
space cannot be paved only with regular tetrahedrons. With mathematical 
rigour it was demonstrated that it was possible only in the case of a 4- 
hypcrsphcrc space [7]. D. Nelson demonstrated that in order to transform 
a hypersphere into a flat space, it was necessary to introduce linear defects 
(like disclination and dislocation) into the structure [8]. Thus, the system 
of linear defects is an integral element of the polytetrahedral structure, and, 
owing to its topological stability, disclination can be considered as a basic 
structural element of the liquid structure. Hence, we can represent the liquid 
structure as a disclination system. 

In order to describe the disclination system, we use gauge theory of 
defects [9]. To write down the Lagrangian of the defect system, let us make 
use of the standard theory of elasticity. The simplest Lagrangian, describing 
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the system with elastic deformations, has the form: 

where Xii^^ ^) is the elastic strain field, A and /i are the Lame constants, po 
is the mass density (which is considered to be constant for simplicity), and 
Uij are the relative deformation components: 

Uab = Cab - Sab = daX'dbX' - Sab 

(Greek letters a, (3, ... arc used to denote the space components set {a} = 
{1, 2, 3}, and Roman letters a, b, ... are used to denote the full index set, 
including the time component {a} = {1, 2, 3, 4}). 

According to the gauge theory of dislocations and disclinations, the plas- 
tic deformation of the matter structure can be considered as breakdown of 
the homogeneity of the rotation and translation (SO(3)l>T(3)) groups action. 
In order to take into account these homogeneity breakdowns, the compen- 
sating fields are introduced into the Lagrangian {A^ and ifl), and transition 
from ordinary to covariant derivatives is carried out: 

daX''-B:, = daX' + ll.,A^x' + ^, 

where Yaj three generating matrixes of the semisimple group S0(3). 
After that the Lagrangian Lq is substituted by the new Lagrangian: 

L = Lo + siLi + S2L2, 

where si and S2 are free parameters of the theory, and the first term describes 
elastic properties of the matter: 

1 j j 1 

Lq = -P0-B4-B4 - - [XEaaEfjfj + 2^Eaj3Eaj3\ 

(where Eab = BiaBu, — Sab is the strain tensor). The second term, 
describes the dislocations, the following denotation is used here: 

Dl, = da^l-db^l+faj - A^^a + F^bX') (F^b = daA^-dbA^+C^^X'^l), 
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in general case Cg^ are constants of the structure of the semisimple group 
S0(3) (C^^ = £^^). The third term, 

S2L2 = -^S2C„/3i^„"65 (5"^ = -<^"^ = 1/0, (1) 

describes dischntions. The Yang- Mills fields, A'^, and describe disclina- 
tions and dislocations accordingly. 

Usually it is supposed that the disclinations are most important struc- 
tural elements in the description of amorphous matter, since they have the 
largest energy density S2 ^ si, and determine two fundamental properties 
of glasses: the absence of long-range ordering and the resistance to crys- 
tallization. It is natural to suppose that disclinations are also fundamental 
structure elements in liquids, therefore hereinafter we will focus only on the 
case of breakdown of homogeneity of the rotation semi-group (monoid). In 
other words, we will neglect the c/p^ field contribution to the action and re- 
strict ourselves to the consideration of the theory with purely disclination 
Lagrangian, 

L = Lq — S2L2. 

In the case of linear disclination the difficulties, concerning nonlinear- 
ity of the Yang-Mills fields, can be avoided. The field of the disclination 
line element differs from the point disclination field, since in this case the 
rotational displacements, corresponding to the gauge transforms, can be per- 
formed only around the tangent to the defect line vector, directed along this 
element (wedge disclination), and, consequently, the gauge group reduces 
to S0(2). Therefore, the theory becomes much simpler since the rotation 
group SO (2) is Abclian instead of the non-Abelian group SO (3) (In work 
[10] this fact was used in the case of infinite rectilinear defects). In order to 
find the field of linear disclination, let us introduce the tensorial fields 
and P^, which take the following form: 

then the energy functional of the "free" disclination field (1) can take the 
form of: 

L2 = \l {d'x)Ft,F-^ = \j {<fx) m - G%m - Gt)] + j {d'x)P^Gl 
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Minimizing this part of the Lagrangian we can get the quasi-stationary ex- 
pression for the potential of the gauge field of the wedge disclination frag- 



menton dla- 



^4 ~ + 2Ci^ dla, ^fc — ^akj^j ^"3 + dla, 

where C and Ci are arbitrary constants. It should be noted that in this 
expression there is no summing with respect to a, and Ci = because of 
A^(oo) = 0. 

Let us consider the two basic examples: 1) First we find the field of 
the infinite rectilinear disclination, directed along the axis ?'„. For that we 
introduce the normal to the axis ia vector Rj, which determines the point 
location with respect to disclination. Then the expression of the field of the 
rectilinear disclination, directed along the axis z, can be written in the form 
of: 



Aj^ — C £zij / Q dlz — C £aij , » q /9 — 2C £zij o • 

J r J ( r>2 , 1 2\'^i'^ R 



Due to the condition $7 = 2t^u = § Afdli, where Q is the Frank vector, and 
u is the Frank index, we get the expression which agrees with the result of 
[10]. 2) The field of the element of the round (with the radius equal to a) 
disclination loop has the form of: 

A^--s -^dl -—n-dl 
z r ^ r 

where rii is the unit vector normal to the loop's plane. 

If we consider only the low-angle disclinations, v ^ 0.2, which arc typical 
for the liquid-like structure, the stress field around the disclination line has 
the form of: 

Uai = llEai = /U {BajBji - 5ai) = /U [SaalXlAf + EailXiA^ + diUa + daUi] . 

Hence, the volumetric density of elastic interaction energy of disclnations, 
described by the fields A and A', has the form of: 

F = ^(JaiEai ~ f^EaalXlAfe^ajx'jA'y + O(a^). 

One can see that form of this interaction is similar to the form of electrody- 
namic interaction of currents. For example, the elastic interaction energy of 
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two round loops F and T' with the radius a, being remote from r, has the 
form: 



r r' 

This expression agrees with the expression of the interaction energy of cur- 
rent loops on condition that n\\n'. 

3 Description of glass transition 

The basic idea of our approach was first suggested by N. Rivier [11]. Accord- 
ing to his theory the glass transition can be considered as phase transition 
in the system of topological defects. In this theory there is no principal 
distinction between the spin or structure-disordered systems. Both in the 
first and in the second cases the problem of description of the system is 
reduced to the description of the system of interacting disclinations. We 
would like to make use of this analogy and reduce this description to that 
of the disordered spin systems as much as possible, and apply the methods, 
known in the spin glass physics, to our system. 

Computer simulation results testify that the disclination system is a tan- 
gled net of interacting disclinations. To describe this system, let us make 
use of the approach, which is well-known in classical electrodynamics [12]: 
we associate the system of linear disclinations, included in some small vol- 
ume, such that the numbers of positive and negative disclinations entering 
this volume are equal, with their general topological moment (the analogue 
to the magnetic moment of the electric currents system). Thus, the total 
elastic interaction energy between all the elements of the disclination net- 
work can be represented as the interaction energy of the system of the local 
topological moments. Since the disclination loop is the simplest structure 
element, having a topological moment, to simplify the model, let us imagine 
the disclination network as the system of randomly located and randomly 
directed disclination loops. 

Let us determine the topology moment vector Si in our system as mag- 
netic moment is determined in electrodynamics. In this case the Hamilto- 
nian, that describes the disclination loops system, can be represented in the 
form of Heisenberg model Hamiltonian: 
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with the inter-spin couphng (from (2)) 

cos{lSiSj] 



J, 



Ik 



,4„6 



/ITT a ViVj- 



where I and k are the indexes of the space coordinates of the vector fij. The 
module and sign of the dischnation loops coupling depend both on their 
situation and their mutual orientation. This system is frustrated. Indeed, 
from the last expression one can see that the moment coupling is an alter- 
nating quantity, hence, with equal possibility these moments can interact in 
either ferromagnetic- like or antiferromagnetic-like manner {{Jij) = 0). This 
character of coupling is typical for spin glasses [13]. Therefore, we attempt 
to describe the considered system, using the modification of the Edwards- 
Anderson model with large but finite-range interaction, which was analyzed 
in [13, 14]. According to this analysis, glass transition in such systems is 
the so-called "hierarchical" phase transition, which is hierarchy of successive 
transitions that closes at the glass transition point Tg. Using the results of 
[14] and the expression of defect coupling energy obtained above one can 
estimate the transition points Tj for every i-th step of this series, as well as 
the correlation lengths corresponding to these points: can be estimated 
from: 



In 



Ti-Tg 




P 


[ ^9 \ 




[ g-1 



iq^+' - 1) 



(3) 



where q = 32/15, p = 4/15 [14], the coordination number, Zq, can be 
evaluated from the coupling radius R: 

where c is disclination density, and R can by determined from the inter- 
disclination interaction energy and thermal fluctuations energy parity con- 
dition: 



E{A) 



yua^vrV^ 3kT 



(4) 



2i?3 2 ' 

where a is a typical interatomic distance. If we express R from (4) and 
substitute this expression in (3) we find that the freezing process begins at 
the temperature, which can be estimated from solving the equation 



1 + ^ 



1 + 



4//C7r'ri*'//2 



(5) 
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Subsequent cooling down results in formation of the clusters which finally 
have shaped at the temperature: 



0.84 




The correlation length characterizing the average size of these clusters as 
well as their average life time can be determined as: 

r Ri tqZq'^^^. (6) 

For quantitative estimation of these parameters let us assume their typical 
values: o = 1 • lO'^ m; k = 1.38 • lO'^^ J/K, Tg = 200 K; tq « 10-^° s; 
the estimation of the disclination density follows from the analysis of the 
experimental data (typical interdisclination distance is ~ 5a): c = 1.4 • 10^^ 
m~^; and the shear modulus of liquids is /U ~ 10^^ N/m^, according to [15]. 
As a result, after substitution of these parameters in (5) and (6), we get 
the following estimations: Tq - Tg ^ 100 K; Ti - w 20 K; ^ w 300 nm; 
T ~ 0.4 • 10~^ s. Of course, these estimations are rather proximate yet. But 
we hope, that they can be improved in the future. 



4 Conclusions 

The structure frustrations and complicated cooperative character of atoms 
moving in liquids lead to the necessity of using of the state-of-the-art meth- 
ods of the statistical physics to understand the nature of its properties. In 
this paper we have tried to combine the gauge fields theory and Anderson 
localization theory to describe the structure formation processes in super- 
cooled liquids. 

The given above estimations allow us to conclude, that the process of 
the freedom degree freezing begins at the temperatures that are significantly 
above the glass transition temperature, Tg, and in the temperature range 
of (T — Tg)~100 K in the structure of the liquid the correlations with the 
size ~100 nm form. The disclination model allows us to explain the density 
beating by natural way, since the atom density is larger close to a negative 
{u > 0) disclination core and smaller close to a positive {u < 0) one. There- 
fore, we think that these correlations conform to the long-range density 
correlations (Fischer clusters), which are observed in the experiments. 

This study was supported by the RFBR grant (04-03-96020-r2004ural), 
and grant of FSP "Integracia" (B0086/2121). 
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